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Precise measurements of the mass and width of the W boson are carried out in e + e~ collisions 
at Lep, by kinematic reconstruction of the invariant mass distributions of W + W~ — > qc{£Ue 
and W + W~ — > qq'qq' candidate events. The most recent combination of such results from 
Aleph, Delphi, L3 and Opal uses approximately 82% of the final Lep2 integrated luminosity 
and is preliminary. The mass of the W boson so determined is Mw = 80.446 ± 0.026 ± 
0.030 GeV, while the corresponding direct measurement of the W boson width gives Tw = 
2.148±0.071±0.063 GeV. These measurements are subject to sizeable systematic uncertainties 
from the QCD phenomena of Bose-Einstein correlations and colour reconnection. Recent 
substantial progress in both of these areas is reported. 



1 W Mass and Width 

The Lep e + e~ collider at CERN has provided an ideal environment for the study of the prop- 
erties of the gauge bosons of the Standard Model of electroweak interactions. Since 1996, it 
operated at centre-of-mass energies above the W + W~ production threshold (Lep2), allowing 
direct measurements of the W boson mass, M\y. When combined with the direct measurements 
of the top quark mass at the Tevatron, these allow further constraints to be set on the mass of 
the Higgs boson via electroweak radiative corrections. Comparison between the direct measure- 
ments of the mass of the W boson and the value determined indirectly from data recorded at 
^ ~ M% provides an important test of the self-consistency of the Standard Model. The direct 
measurement of Tw further tests the consistency of the Standard Model. 

The measurement of M\y an d Tw - is divided into three stages: selection of W + W~ events, 
event-by-event mass reconstruction, and determination of M\y itself. 



1.1 W+W" Final States 



W + W~ events are divided into three final states. W + W~ — > qq'qq' events comprise 45% of 
the total W + W~ cross-section and are characterised by four energetic jets of hadrons with 
little or no missing energy. Semi-leptonic W + W~ — ► qq'tPe decays comprise 44% of the total 
W + W~ cross-section and are characterised by two distinct hadronic jets, a high- momentum 
lepton and missing momentum due to the prompt neutrino from the leptonic W decay. The 
signature for the W + W~ — > qq'r^ r channel is similar, with the exception that the r lepton is 
identified as an isolated, low-multiplicity jet typically consisting of one or three tracks. The 
W + W~ — ► £ + vg£'~Uj>i channel, with at least two unobserved neutrinos and a relatively low 
branching fraction, has limited Mvy sensitivity and is not discussed herein. 

1.2 Invariant Mass Reconstruction 

The clean environment at Lep allows a complete kinematic reconstruction of the invariant mass 
on an event-by-event basis. Hadrons are grouped together into jets using clustering algorithms 
such as k±. In W + W~ — * qq'tPi events, charged leptons are identified and neutrinos are inferred 
from the missing energy and momentum. The invariant masses of the two W bosons can be 
determined directly from the reconstructed momenta of observed decay products. Experimen- 
tally, the limiting factor in the mass resolution is the uncertainty in the jet energy measurement, 
which is poor in contrast to the measured jet directions. As the centre-of-mass energy is well 
known, the mass resolution can be improved significantly (factor ~2-3) by use of a constrained 
kinematic fit imposing the four constraints of energy and momentum conservation (4-C fit). 
Small additional gains are possible by imposing the additional constraint that the masses of the 
two W bosons are equal in each event (5-C fit), giving a single mass measurement per event. 

For W + W _ — > qq'iVe events, the number of effective constraints is reduced to 2 (1) for 
a 5-C (4-C) fit due to the three missing degrees of freedom corresponding to the unmeasured 
neutrino momentum. For W + W~ — * qq'r^ T events, most of the mass information is given by 
the hadronically decaying W. A frequent assumption in constructing kinematic fits for these 
events is that the true r direction coincides with its observed decay products, while the r energy 
is unknown, removing a further constraint. 

In W + W~ — * qq'qq' events where four jets are reconstructed, there are three possible pair- 
wise combinations. In some analyses, a single preferred combination is selected on the basis of 
information such as kinematic fit probabilities, jet-jet angles or the CC03|] matrix elements. In 
other analyses, all combinations are used with different weights assigned to each, or two combi- 
nations are used with equal weight. In addition, as quarks may radiate energetic gluons leading 
to a distinct five-jet topology, several analyses separate events into four- and five-jet categories 
to be treated separately, leading to an overall improved mass resolution. 

1.3 M\y Determination 

There are three main methods whereby Jl% is determined. The most widely used (Aleph, L3, 
Opal) involves reweighting Monte Carlo events including detector simulation to an arbitrary 
value of M\y using the ratio of 4-fermion or CC03 matrix elements. A likelihood fit to the 
reweighted mass spectra determines the value of J\% which best resembles the data. The fit can 
be either a 1-dimensional fit to the 5-C fit mass, possibly performed in several exclusive regions 
of reconstructed mass error, or a 2-dimensional fit to the 4-C fit masses. There is an implicit 
MC correction, e.g. for effects such as initial state photon radiation (ISR) and detector effects. 

A second method (Delphi, Opal) constructs an event likelihood from the convolution of a 
Breit-Wigner with a radiator function (to account for ISR) and a resolution function (to account 

a Doubly-resonant W + W~ production diagrams, i.e. t-channel v c exchange and s-channel Z°/7 exchange. 



for detector effects). The probability of an event being W + W~ is also included in the likelihood. 
These 'convolution' based methods make more use of the information per event and so might 
be expected to give a more precise Jl% determination. An explicit bias correction has to be 
applied to the fitted Mw, based on events including detector simulation. 

The third method (Opal) fits an analytic function (asymmetric Breit-Wigner) to the re- 
constructed mass spectra. Similarly, an explicit bias correction is made. All methods yield 
comparable precision on the measured Mw- 

1.4 Systematic Effects 
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Table 1: Summary of the systematic uncertainties for the combined Lep fit results. 

The statistical uncertainty on the combined Lep Mw measurements is 26 MeV, therefore a 
clear understanding of systematic uncertainties is vitally important. Systematics may be corre- 
lated between combinations of experiments, channels and years of data taking. A summary of 
the uncertainties, without giving the detailed decomposition into the various correlated sources, 
is given in Table |l|. A consequence of the large systematics from Bose-Einstein Correlations 
(BEC) and Colour Reconnection (CR) is that the W + W~ — > qq'qq' channel has a reduced con- 
tribution to the overall Mw, with a weight of 0.27. In the case where both channels had equal 
weight, the statistical uncertainty on Afw would be reduced by ~ 15%. 

The most significant systematics are those which are correlated between experiments, as 
described below. Photonic radiative corrections (ISR, FSR) have been estimated originally by 
comparing different ISR models, e.g. as in Koralw and Excalibur. More recently, Koralw 
events have been reweighted to correspond to an 0(a 2 ) or 0(a) treatment of ISR using the 
matrix elements calculated inside the model (the default calculation is at C(a 3 )). However, 
these calculations are incomplete at 0(a) (no ISR/FSR interference, or direct 7 radiation from 
Ws) , so ideally a complete calculation as in the Double Pole Approximation of Racoon WW would 
be used. As yet, this only produces weighted events so no realistic estimates are available. 

Hadronisation has been estimated by comparing models, e.g. Jetset and Herwig. Other 
approaches include reweighting key variables in MC to correspond to data and propagating the 
effect to Mw- Delphi also advocate the use of so-called "mixed Lorentz boosted Z°" (MLBZ). 

The relative uncertainty in the Lep beam energy enters directly into the uncertainty in Mw, 
due to the use of kinematic fits which include a constraint to the centre-of-mass energy. Uncer- 
tainties in beam energy are taken directly from the detailed studies of resonant depolarisation, 
NMR probes/flux loop measurements, and the Lep spectrometer project. 

A significant bias to the apparent W mass measured in the W + W~ — * qq'qq' channel could 
arise if the hadronisation of the two W bosons is not independent and correctly modelled. Final 



state interactions such as CR and BEC may cause just this effect and are estimated by using 
phenomenological models. Direct searches for these effects, which may limit the viable set of 
such models, are described below 

2 Colour Reconnection 

In W + W~ — > qq'qq' events, the products of the W decays in general have a significant space-time 
overlap as the separation of their decay vertices is small compared to characteristic hadronic 
distance scales. Colour reconnection refers to a rearrangement of the colour flow between the two 
W bosons. The effects of interactions between the colour singlets during the perturbative phase 
are expected to be small, ~ ( M " s ) 2 Tw. The situation is less clear in the non-perturbative 

*l * colour s 

phase, where phenomenological models are implemented in hadronic Monte Carlos. A higher 
susceptibility to CR (and more Z°/7 — > qq background) is expected when W + and W~ hadroni- 
sation regions overlap, so the spacetime picture of the QCD shower development is important. 

The predicted (barely) observable effects of CR include changes to the charged particle 
multiplicity, momentum distributions and the particle flow relative to the 4-jet topology. The 
aim is to establish whether CR actually exists as well as controlling (or better, calibrating using 
data) the bias on M\y measurements. The basic analysis method consists of comparing fully 
hadronic events with either: models with/without CR; W + W~ — > qq'iUi events; or MLBZs. 

2.1 Charged Particle Multiplicity 
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Table 2: Summary of recent charged particle multiplicity measurements, with combined statistical and systematic 
errors. Aleph measurements are not corrected for detector effects. 

The difference in multiplicities between W + W~ — ► qq'qq' and twice the hadronic component 
of W+W- -► qq'^ events, A(n ch ) = (njjf) - 2(n%f u ), should be zero in the absence of CR 
effects. A compilation of the most recent multiplicity measurements is given in Table The 
effects are expected to be enhanced for softer particles, p < Tw, so observables with implicit 
scale such as ln(l/x p ), px, rapidity, have also been studied with similarly inconclusive results. 
As the systematics are now comparable to the statistical uncertainties, no attempt is made to 
combine results in a trivial way. Delphi and Opal have studied the effects of CR for heavier 
hadron species, such as kaons and protons. Although the numerical effects of CR are larger, the 
significance is reduced due to loss in statistical precision. 



2.2 Interjet Analysis 

These analyses §,!,§, are based on the "string effect" analysis in e + e~ — ► qqg at lower energies. 
All Lep collaborations are now using these analyses, recently developed by L3 and following 
earlier studiesB. Events are selected by requiring, in addition to usual W + W~ — ► qq'qq' criteria: 
four distinct jets (1/34 > 0.01 in k± scheme); the two largest jet-jet angles are 100°-140° (intra- 
W regions); the two smallest jet-jet angles are < 100° (inter-W regions) and have no jets in 
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Figure 1: Particle density in a single di-jet 
plane. 



Figure 2: Particle density for all planes, 
scaled to di-jet opening angles. 



common. This leads to an efficiency ~ 15%, good jet-jet to W association, with "strings" 
that are back-to-back and not crossing one another in 87% of events. In the case of L3 at 
y/s = 189 GeV, 209 such events are selected. Next, particles are projected onto each of the four 
selected di-jet planes, forming the particle density as a function of angle from one of the plane 
defining jets. This is illustrated in Figure [l], where the inter-W regions correspond to the larger 
regions between the jet peaks. As the di-jet angles vary from event to event and plane to plane, 
particle densities are scaled to the di-jet opening angle event-by-event, as show in Figure ||, after 
background subtraction. 
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Figure 3: Ratio of inter-W/intra-W particle flow (a) 
higher efficiency analysis, weighted by In — , (b) L3-like 

selection, seeB for details. 



To quantify the effects, the ratio of the inter- 
W/intra-W regions (or its inverse) is formed, 
as shown in Figure [3|, and integrated in the 
central region away from the jet peaks. An 
alternative is to consider the ratio of the in- 
tegrals of inter-W and intra-W particle den- 
sities. L3 estimate their sensitivity to the 
SK I CR model!, is 3.2 (0.5) a total error 
for 100 (32) % reconnected events, using | 
of their Lep2 sample. Opal have a variant 
on the analysis using a likelihood based se- 
lection to associate di-jets with Ws, a 4-C 
kinematic fit to define jet axes, and no dou- 
ble counting of particles, which predicts a 
slight improvement in sensitivity and favours 
~ 65% of events reconnected in the SK I 
model. However, their emulation of the L3 
analysis prefers a no-CR scenario. 



2.3 Towards Mass Biases 



In the SK I model, the reconnection probability is a free parameter and therefore can be adjusted 
to provide the best agreement of data. L3 estimate their data prefer ~ 40% CR in this model, 
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Figure 4: x between data and SK CR model as a function of reconnection probability. 



and are 1.7a separated from no-CR scenario, as shown in Figure |3|. A similar study in Aleph 
suggests data prefer this model dependent probability to be ~ 15%, also shown in Figure |||. 
Having estimated the preferred CR probability, it is trivial to estimate the corresponding bias 
on Myy in the context of a single model. Comparison of such results between collaborations is 
complicated by differences in the models they use, such as hadronisation tuning. Using common 
samples of simulated events which are processed by each collaboration's detector simulation is 
in progress and will lead to a significant improvement in understanding in this area. 



3 Bose-Einstein Correlations 



Bose-Einstein correlation leads to the enhanced production of identical boson pairs, such as 
7r + 7r + or 7r~7r~, at small 4-momentum difference, Q 2 = — (pi — P2) 2 - This phenomena is firmly 
established in various environments, in Z°/7 — > qq at LepI and between the particles of a 
single W boson at Lep2, among others. Traditionally, it is studied using a 2-particle correlation 
function: i?(pi,p2) = p2i.P1, P2)/po(Pi) P2), where P2 and po are 2-particle densities with and 
without BEC, respectively. Of particular interest in the W + W~ case is whether or not there is 
any additional correlation introduced between the decay products of the two Ws, which could 
potentially bias M\y measurements. 

One serious problem in this area is the construction of the reference sample, pq, and there 
are three frequently used methods. The first takes unlike-sign particle species, such as ir + ir~ , 
forms the ratio of like-sign/unlike-sign, and then takes the ratio of this quantity relative to 
Monte Carlo to reduce the impact of having resonances in po but not in p2 ("double ratio"). The 
second consists of taking po to be a like-sign MC sample, although this is subject to deficiencies in 
modelling. The third involves mixing pairs of data events, such as Z°/7 — > qq or the hadronically 
decaying W in W + W~ —* qq'iUi. Another serious problem is the ignorance of non-perturbative 
QCD amplitudes, forcing analyses to resort to MC models, which in turn suffer from being 
probabilistic in nature, with the models of BEC implemented in Monte Carlo simulations. 

A common parametrisation of the correlation function is R(Q) ~ 1 + A exp(— r 2 Q 2 ), where 
A and r represent the source strength and size, respectively. Results from each of AlephS, L3 
and Opalq are summarised below. No Delphi results were presented (at their request), as 
they were undergoing substantial revision. Results which differed qualitatively from their earlier 
analyses were presented at Moriond QCD 2001. 

In the Opal analysis, the two-particle correlation function is formed using the double ratio 
of like-sign/unlike-sign, data/no-BEC MC. Three event classes are identified, W + W~ — > qq'qq' 



W + W~ — > c\q'£v£ and Z°/7 — ► gg, each consisting of the linear sum of pure contributions 
weighted by a probability derived from MC. The classes are C DIFF (inter-W BEC), C SAME 
(intra-W BEC), C z * (non-radiative qq. A simultaneous fit is performed to extract \ DIFF for 
various source size hypotheses. Examples are: R DIFF = R SAME = R z * j yielding \ DIFF = 
-0.14±0.36, and completely independent Rs, giving \ DIFF = 2.9+1.7 and \ SAMF = 0.62±0.10. 
Intra-W BEC is established, but the analysis is unable to ascertain whether inter-W BEC exist. 

In the Aleph analysis, again the double ratio is used, and the Z°/7 — ► qq background, in 
which BEC are known to be present, is modelled using the model BE%. It is concluded that 
inter-W BEC is disfavoured by data with a significance of 2.2<t, while the data are compatible 
with intra-W BEC. An event mixing analysis is also used, which qualitatively disfavours inter-W 
BEC, but is as yet incomplete in that a full study of systematic errors has not been performed. 

The L3 analysis uses the method of Chekanov et alu, which sets out a very robust framework 
to test for the presence of inter-W BEC. If the W + and W~ decays are uncorrelated, then: 
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where the first two terms on the r.h.s. are estimated from individual W + W~ — > qq'fl^ events, 
assuming they are the same in W + and W~ decays, while the rightmost term is formed by 
mixing pairs of W + W~ — ► qq'-fi^ events to give "fake" events in which there can be no true 
inter-W BEC. All background samples include BEC using the BE32 model. 
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The variable D is defined as the ratio of the 
left hand side to the right hand side of Equa- 
tion [l]. In the absence of inter-W correlations, 
such as from BEC, and of bias introduced in 
the event mixing, D = 1. To eliminate this 
potential residual experimental bias, the vari- 
able D' is constructed as ratio of D in data to 
that in MC (having only intra-W BEC). Fig- 
ures |B|(a)-(b) and |5|(c)-(d) show the variables D 
for like-sign and unlike-sign data, and similarly 
for D' . The data clearly favour the no inter-W 
BEC model. By fitting the phenomenological 
correlation function to these data, L3 obtain 
A = 0.013 ± 0.018 ± 0.015, where A = cor- 
responds to no inter-W BEC. In contrast, the 
BE 32 model gives A = 0.126±0.006(stat.), thus 
the data disfavour inter-W BEC by 4.7cr. 



Figure 5: Inter-W BEC observables, see text for details. 



4 Combined Lep Results and Summary 

Based on an analysis of 2.8 fb _1 (~ 82% of the entire Lep2 data sample, 100% analysed by Aleph 
and L3), the following preliminary measurements!, summarised in Figures ^ and [/], are made: 
M w = 80.446 ± 0.026 ± 0.030 GeV, T w = 2.148 ± 0.071 ± 0.063 GeV, A A / W (qq'qq' - qq'^i) = 
+18 ± 46 MeV. Further improvements in systematics are anticipated. Colour reconnection 
analyses, such as interjet multiplicity, from all collaborations will be combined. The overall 
conclusion regarding inter-W Bose-Einstein correlations among the Lep collaborations is finally 
becoming more consistent: there is increasingly strong evidence that they do not exist. 
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Figure 6: Combined Lep Mw and Tw, all channels. 
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Figure 7: Comparison of Lep Mw with other measurements and Standard Model predictions 
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